Can one hear the density of a drum? Weyl's law for inhomogeneous media 
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We generalize Weyl's law to inhomogeneous bodies in d dimensions. Using a perturbation scheme 
recently obtained by us in Ref. we have derived an explicit formula, which describes the asymp- 
totic behavior of the eigenvalues of the negative laplacian on a closed d- dimensional cubic domain, 
either with Dirichlet or Neumann boundary conditions. For homogeneous bodies, the leading term in 
our formula reduces to the standard expression for Weyl's law. We have also used Weyl's conjecture 
to obtain a non-perturbative extension of our formula and we have compared our analytical results 
with the precise numerical results obtained using the Conformal Collocation Method of Refs. [3, @]. 



According to Weyl's law the largest frequencies of the 
sound of a uniform drum (membrane) are primarily de- 
termined by the area of the drum and not by its shape. 
The same result also applies to the electromagnetic field 
inside a waveguide, or to the cigenmodes of a particle 
in a quantum billiard, or more in general to the normal 
modes of the Laplacian operator in a closed domain in 
d dimensions, with Dirichlet or von Neumann boundary 
conditions. 

Weyl conjectured that the number of modes in a two 
dimensional drum of area A and perimeter L, of energy 
lower than a given energy E is 

N(E) = —E^—^JE+... 

V ' 47T 47T 

where the — and + signs hold for Dirichlet and Neumann 
be respectively. This formula is known as Weyl conjec- 
ture and, if true, it implies 
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which describes quite well the spectrum of a homoge- 
neous drum. This conjecture has also been extended to 
d dimension. Ref. which contains a detailed account 
of Weyl's law, provides a historical perspective. 

In this paper we want to ask ourselves a question sim- 
ilar to the one that Weyl posed himself: how will an 
inhomogeneous drum, with density varying from point 
to point, sound? or, more in general, what is the asymp- 
totic behavior of the eigenvalues of the inhomogeneous 
Helmholtz equation on a d-dimcnsional domain, with a 
spatially depending density, and with Dirichlet or Neu- 
mann boundary conditions? And, finally, is it possible to 
justify Weyl's conjecture using perturbation theory? 

To the best of our knowledge this problem has not been 
formulated before Q although its solution could find an 
incredible number of applications: for example to study 



1 In Ref. 3 the inhomogeneous media problem is briefly mentioned 
in the section concerning "open problems". 



of the acoustics of an inhomogeneous medium, or the 
properties of an electromagnetic wave-guide with a vary- 
ing dielectric coefficient, or even to study the propagation 
of a seismic wave traveling in the interior of the earth. An 
asymptotic Weyl-type formula for inhomogeneous bodies 
would also be helpful to extract informations on the den- 
sity of the material composing the body, thus allowing to 
rephrase Kac's question [B| "can one hear the shape of a 
drum?" into "can one hear the density of a drum?" 

Although these are difficult questions, we will show in 
this paper that it is possible to answer them, obtaining 
accurate estimates. We hope that the formalism devised 
here will pave the way to more refined calculations and 
also provide an alternative way to look at the old and 
important problem considered by Weyl. 

Let us describe our approach. We focus for the moment 
being on two dimensional homogeneous membranes of ar- 
bitrary shape. In a recent paper, Ref. we have devised 
a perturbative approach to the calculation of the normal 
modes of these membranes (or quantum billiards) , that 
involves a conformal map which sends the original shape 
into a reference shape (square, circle, etc.). 

To be more explicit, w = f(z) is a conformal transfor- 
mation that maps a region f2, which for example could 
be chosen to be a circle or a rectangle, into the region T>, 
representing the shape of the drum. Under such transfor- 
mation the homogeneous Helmholtz equation transforms 
to an inhomogeneous Helmholtz equation: 



Aip(x, y) + EE(x, y)ip(x, y) = , 



where 



E = 



(2) 



(3) 



We will refer to E as the conformal density. Notice that 
E could also be interpreted from the start as a physical 
density, although the reverse is not necessarily true, as 
we shall soon see. 

If we have in mind shapes which are obtained by small 
deformations of the reference shape f2, we may express 
E = 1 + <t, where a is the perturbation density generated 
by the mapping. In Ref. [l[ we have obtained an explicit 
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form for the perturbativc corrections to the energy up to 
third order in a: 
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Clearly e„ and |n) are the exact eigenvalues and eigen- 
states on f2 (n is the set of quantum numbers which define 
the state). We have defined u n k = e« — (k- 

As we have observed in Ref. [l| the first terms ap- 
pearing in each of the equations above correspond to the 
terms of a geometric series: 

e„ (1 - (n\a\n) + (n\a\n) 2 - (n\a\n) 3 + .. .) , 



which can be resummed as 



1 + (nan) 
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(8) 



It is important to realize that, although the interpre- 
tation of S as a conformal density is limited to two di- 
mensional problems, the perturbative scheme of Ref. [l[ 
is general and it can be applied to a problem in d dimen- 
sions, where £ is a physical density. In what follows we 
will therefore work in d dimensions, assuming to be a 
d-dimcnsional cube of side 2L centered in the origin and 
£ a physical density inside fid- 

The eigenfunctions of fid corresponding to Dirichlet 
and von Neumann boundary conditions are obtained with 
the direct product of the functions on each orthogonal 
direction. We have 



*£)( x ) = -L S m[^(x + L) 



for the Dirichlet modes (n x = 1,2, 



and 



= cos 

™ x y/Z 
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for the von Neumann modes (this expression holds for 



1,2, 



for n. 



0, * 



= 1/V2L). 

Let us now go back to eq. JS]): we first concentrate 
on the matrix elements and look for a suitable 

approximation for n — > oo. As we have already men- 
tioned here n stands for the full set of quantum numbers 
specifying the states, i.e n = (n Xl , . . . ,n Xd ). 



Therefore 

(n|S|n)= / ¥£f)\x 1 )...*£p 2 (x d )X(x 1 ,...,x d ) 
For n Xi — ► oo one may approximate the highly oscil- 

f D TV) 

latory functions with their average value, $„,.' (xi) « 
1/2L, and therefore write 

(n|S|n) w S(a;i,.-.,^) ■ 

We now would like to relate the energy of a state in 
fid, with quantum numbers n Xl , . . . , n x 
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to the number of states with equal or lower energy N(e). 

As shown in Ref. [f| the number of states with Dirichlet 
be of energy less than E na . na . may be expressed as 
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For states obeying von Neumann be we have 

N {N) {E) ^ N {D) (E) 

2 {2LY-H^ /2 _ 1/2 
(4^) d / 2 r(d/2 + 1/2) ' 1 J 

where the second term counts the states with one of the 
quantum numbers vanishing. 

Substituting these results in eq. ([8]) we obtain a relation 
for the energy as a function of N: 
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which is a generalization of Weyl's law for a d-cube filled 
with density E(xi, . . . , Xd)- The + and — signs hold for 
Dirichlet and Neumann be respectively. 

It is interesting to check some particular limits of this 
expression. For example, for £ = 1, one should recover 
the standard form of Weyl's law for an homogeneous d- 
cube: in this case the first term of eq. (fT3)) reduces to the 
correct expression 
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where Vd = (2L) d is the volume of the d-cube. 

Let us now focus on the case d = 2. In this case we 
have that 
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which should be compared with eq. {T]). Assuming that 
E is now the conformal density previously discussed we 



may state Weyl's conjecture, eq. (p}, in the form 
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J n ^(x,y)dxdy J Q T,(x,y)dxdy y J n ^(x,y)dxdy 

I 



(16) 




FIG. 1. (color online) Cardioid drum with density p(u,v) 
1/(1 + 4(u 2 + v 2 )). 
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FIG. 2. (color online) 

log 10 |l- [f n X(x,y)dxdy] (4f> 



The quantity Hjv = 
for the 



+ E ( N N ^)/8nN 



cardioid drum of Fig. [T] The energies ' and have 
been calculated using the Conformal Collocation method of 
Ref. 0] with a grid of 99 points on each axis. 



after noticing that J n T,(x,y)dxdy and that 
E 1/,2 (x, y)ds are the area and perimeter of the 
membrane. The reader may observe that eq. Q16p 
reduces to our eq. (fTS*)) in the " perturbative" regime, 
\<t\ -C 1. As a matter of fact, in this limit one has that 

J gn E 1 / 2 ^, y)ds/ J n E(x, y)dxdy m 4 (the perimeter to 

square root of the area ratio for a square), from which 
eq. (fl~5j) follows. 

We will now use both eqns. |T]) and (fT5)) to study two 
dimensional drums of arbitrary shape and density. The 
transformed Helmholtz equation now reads 

Atp(x, y) + EE(x, y)p(u, v) ip(x, y) = , 

where p(u, v) is the physical density of the drum and 
(u, v) = (Ref(z),Imf(z)). Therefore eqns. ([T]) and (fl~5|) 
must now be used substituting E(ai, y) with E(x,y) = 
E(x,y)p(u,v). 

We have tested eq. ((T5j) and (p~6|) on a cardioid drum 



with density p(u, v) = l/(l + 4(u 2 + u 2 )), as seen in FigfJ] 



In this case we have that 



E 1/2 (a;,y) ds « 3.00112 



A 



E(x,y) dxdy w 1.21205 



(17) 
(18) 



Notice that L and A do not have a geometric interpre- 
tation of perimeter and area of a drum since the ratio 
L/A = 2.48 is smaller than the corresponding ratio be- 
tween circumference and area of a circle with area A, 
2/VttA = 3.22. Therefore E(x, y) cannot be obtained 
from a conformal transformation. 

An independent confirmation of this observation comes 
from the Payne-Polya- Weinberger conjecture 0,1 (later 
proved by Ashbaugh and Benguria in Ref. Q ) , according 
to which the ratio between the first two eigenvalues of 
the Dirichlet laplacian is maximal for the circle, i.e. 
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2.539 



(19) 



where jo.i and j\ t i arc the first positive zeroes of the 
Bessel functions Jq(x) and J\(x). 
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FIG. 3. (color online) The quantity 6 N = (E^ 1 - >)/2 for 
the cardioid drum of Fig. [1] The fluctuating (blue) line is the 
numerical result obtained using the Conformal Collocation 
method of Ref. @] with a grid of 99 points on each axis; the 
dashed (red) line is the second term in eq. (|16|l while the 
dotted (orange) line is the result obtained eq. Q13JI. 



made by Gordon, Webb and Wolpert in 1992, Ref. [lj|, 
of two different homogeneous drums which are isospec- 
tral, it is now known that the answer to Kac's question 
is negative. In the present case, it make sense to wonder 
if there are truly genuinely inhomogeneous drums (i.e. 
inhomogeneous drums which cannot be reduced to a ho- 
mogeneous one by means of a conformal transformation), 
which arc isospcctral although corresponding to different 
(inequivalcnt) domains and/or densities. 

We conclude by observing that an improvement of the 
analytical formulas contained in this paper necessarily 
involves taking into account the terms in the perturba- 
tion expansion which have been neglected here. As these 
terms involve sums over internal states, we expect them 
to be more difficult to evaluate; moreover, in this case the 
possible degeneration of the levels must also be taken into 
account, as discussed in Ref. Most likely the calcula- 
tion of the corrections to the asymptotic formula of the 
present paper requires some averaging procedure over the 
population of highly excited states. We plan to consider 
this issue in future works. 



Using the Conformal Collocation Method (CCM) with 
a grid of 99 points in each direction we have found 
for the inhomogeneous Hclmholtz equation with density 
p(u,v) = 1/(1 + 4(w 2 + v 2 )) the values 

E{ D) = 10.6769 , E { 2 D) = 29.7008 , (20) 

corresponding to a ratio r = E^ /e[ D ^ m 2.78. Since 
this ratio violates the Theorem proved by Ashbaugh and 
Bcnguria, S cannot be interpreted as a conformal density. 
In other words, it is not possible to build an homogeneous 
drum of appropriate shape so that it sounds precisely as 
this inhomogeneous drum. 

As we see in FigJ5] and [3J eq. (TT6")) describes quite 
precisely the behavior of the energies of this inhomo- 
geneous drum. In particular from Fig. [2] we learn that 
(E^ + E^)/2 essentially behaves as the first term in 
cq. p6[) . i.e. Wcyl's law for inhomogeneous drums. 

For bodies of dimension d > 3, our formulas only apply 
at present to cubical shapes of arbitrary density, since it 
is not clear to the author if a generalization of a conformal 
transformation to higher dimensions exists. However, if 
this extension exists, there is an angle-preserving map 
which relates an arbitrary d-dimcnsional region to a d- 
cube. In this case the formula obtained in the present 
paper straightforwardly applies. 

The results obtained here also allow to give a sense to 
the question contained in the title: after the discovery 
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